Multibreathers in Klein-Gordon chains with interactions beyond nearest neighbors 
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, We study the existence and stability of multibreathers in Klein- Gordon chains with interactions 

that are not restricted to nearest neighbors. We provide a general framework where such long 
range effects can be taken into consideration for arbitrarily varying (as a function of the node 
distance) linear couplings between arbitrary sets of neighbors in the chain. By examining special 
case examples such as three-site breathers with next-nearest- neighbors, we find crucial modifications 
' to the nearest-neighbor picture of one-dimensional oscillators being excited either in- or anti-phase. 

' Configurations with nontrivial phase profiles emerge from or collide with the ones with standard (0 

or tt) phase difference profiles, through supercritical or subcritical bifurcations respectively. Similar 
bifurcations emerge when examining four-site breathers with either next-nearest-neighbor or even 
interactions with the three-nearest one-dimensional neighbors. The latter setting can be thought of 
as a prototype for the two-dimensional building block, namely a square of lattice nodes, which is 
also examined. Our analytical predictions are found to be in very good agreement with numerical 
' results. 

> 

Q\ • I. INTRODUCTION 

' The initial numerical inception of anharmonic modes consisting of a few excited sites in nonlinear lattices |31l . 
■"sj" , [s^], and their subsequent placement on a rigorous existence basis (under rather generically satisfied non-resonance 
■ conditions) in [2^ has triggered a huge growth of interest in the theme of the so-called discrete breathers. These 
. are exponentially localized in space, periodic in time states which have subsequently been theoretically/numerically 
predicted and experimentally verified to arise in a very diverse host of applications. These include (but are not 
limited to) DNA double-strand dynamics in biophysics |33| . coupled waveguide arrays and photorefractive crystals in 
nonlinear optics 0, [l9l . [25j , breathing oscillations in micromechanical cantilever arrays [34] , Bose- Einstein condensates 
^ in optical lattices in atomic physics |30| , and granular crystals [sH] . The interest in this theme has also been mirrored 
, in a wide array of reviews on methods of identifying and analyzing such intrinsically localized modes [1, [T^, [l^, [2^ . 
More recently, the stability of the discrete breather configurations, especially in the case of the excitation of multiple 
sites has been of particular interest. One approach to obtaining relevant results consists of the so-called Aubry band 
theory used e.g. in [1,0. This led to the conclusion that for soft nonlinear potentials multi-breathers with any 
subset of adjacent sites being excited in-phase are unstable, while ones with all adjacent sites in anti-phase can be 
stable in the vicinity of the so-called anti-continuum limit of uncoupled anharmonic oscillators. A complementary 
theory that yields insights on both the existence and the stability of multibreathers has been pioneered by MacKay 
and collaborators; see e.g., 0, [13, [l^- This is the so-called effective Hamiltonian method which is identified by 
averaging over the period of the unperturbed solution and developing the proper action-angle variables. The extrema 
of the resulting effective Hamiltonian determine the relative phases of adjacent excited sites in the multi-site breather 
solution, while the relevant Hessian is intimately connected to the Floquet multipliers of the associated periodic 
orbit. Using this methodology, the work of 22] retrieved as well as refined the results of ^ for arbitrary phase 
relations between the excited oscillators of such multi-breather configurations. The equivalence between these two 
basic methods and their conclusions was recently established in [ll|. We should also note in passing that similar 
results have been acquired also in configurations where there are "holes" between the excited breather sites [3^ . 
through higher order perturbation theory generalizing the above conclusions to the cases with one-site holes. On the 
other hand, the existence and stability of single/multi-site breathers have been studied in diatomic FPU lattices. The 
work of 38| was based on a discrete Sturm theorem which necessitated (for the separation of the space n and time t 
variables) a potential which was at least purely quartic. In the realm of lattices with longer than the nearest-neighbor 
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interactions a variety of issues have been considered such as, e.g. in |14| . the existence and bifurcation of quasi periodic 
travehng waves in nonlocal lattices with polynomial type potentials. 

In the present work, we consider the generalization of the above settings, which are principally concerned with 
the interaction between nearest neighbors, to the case with longer range neighbor interactions for Klein-Gordon 
chains. Upon revisiting the nearest neighbor case and presenting the effective Hamiltonian formalism (of MacKay 
and collaborators) there (section II) for existence and stability of multibreathers, in section III, we generalize this 
formalism to the case of an arbitrary number of neighbors (denoted by r) interacting with each other. By specializing 
to the case of nearest and next-nearest neighbor interactions (and three-site breathers) as our first case example of 
the application of the results in section IV, we already infer the fundamental modifications to the standard picture 
that ensue due to interactions beyond nearest neighbors. These include configurations that have non-standard relative 
phases between adjacent oscillators, a feature which is absent in the nearest-neighbor interaction case [20j and also 
symmetry breaking bifurcations that arise due to the "collision" of branches of solutions with such non-trivial phase 
relations, with more standard ones with relative phases of or tt between adjacent oscillators. The generic nature 
of these conclusions is confirmed by considering the case examples of four-site breathers with next-nearest-neighbor 
interactions in section V and such breathers with interaction ranges of r = 3 in section VI. The latter setting is very 
close to genuinely two-dimensional setting in a square lattice plaquette, which constitutes our final example in section 
VII. We close our presentation by some remarks on the parallels of our results with the simpler case of the discrete 
nonlinear Schrodinger lattices [l3| (section VIII) , which has been examined earlier in 0, [3 1 '^^11 as a summary of 
our conclusions and some future directions (section IX). 



II. BACKGROUND: THE CLASSICAL KLEIN-GORDON CHAIN 




FIG. 1: The classical nearest-neighbor Klein-Gordon chain 

The Hamiltonian of a Klein-Gordon chain with nearest neighbor interactions is the following 

1 



H = Ho+eHi= J2 

i— — o 

which leads to the equations of motion 



+ V{x,) 



(1) 



Xi = -V'{xi) + £(xj_i - 2xj + Xj+i). 

It is well known that this system supports discrete breather, as well as, multibreather solutions. As indicated above, 
there are several papers dealing with the existence and stability of these motions; see e.g. J,, ^, QA^ i32|, i36|] . 



A. Persistence of mutibreathers 



In the anti-continuum limit e = 0, we consider all the oscillators of the chain at rest except for n + \ "central" 
ones which move in periodic orbits of frequency uj. As indicated in 32], it is possible to generalize considerations to 
the case where not all of these oscillators are adjacent to each other, however, we will not concern ourselves with this 
additional complication herein. The time-periodic and space-localized motion of our excited oscillators will persist 
for e ^ to provide multibreathers of the same frequency uj, if the phase difference between the succesive central 
oscillators satisfies specific conditions. In [^l it was shown that multibreathers correspond to critical points of 
which in first order of approximation is given by — Ho{Ii) + e{Hi){<j)i, li) The variables <j)i = Wi+i — Wi 

denote the n phase differences of the n + 1 successive central oscillators, while li are given by — X]j=i ■^j^ where 
( Ji, Wi) are the action-angle variables of each uncoupled oscillator. 

The average value of the coupling part of the Hamiltonian 

(i?i)(0„/,) - ^ <L Hi{wQ,dp,,U)dt 
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is calculated along the orbits in the anti-continuum limit e = 0. 

This yields the conclusion that the persistence conditions for the existence of n + 1-site multibreathers are 



d{H 



il=0, i = (2) 



Note that the persistence conditions are the same for every lattice case where the Hamiltonian can be written in the 
for H = Ho + sHi with ^ 0. For a more detailed description of the above procedure one can also see [22| . 

The motion of the central oscillators for e = can be described by 

oo 

Xi{wi) = ^ Am{Ji)cos{mwi). (3) 

Since the action Ji remains constant along an orbit in the anticontinuum limit, Xi depends only on Wi. So, the average 
value of Hi becomes (|22j] appendix A) 

oo n 
m—1 s—1 

and the persistence conditions ^ become in the case of Klein-Gordon chains with nearest neighbor interactions, 

= ^ Af (<?!>) EE V mAl, sin(TO0,) = Q, i = l...n. (4) 

m— 1 

The function M{(j)) possesses the obvious solutions (pi = 0, tt, while it has no others, as it is shown in [20j . 



B. Stability of multibreathers 



The spectral stability of the above mentioned multibreather solutions or, equivalently, the linear stability of the 
corresponding periodic orbits is determined through its characteristics exponents Gi. These exponents are connected 
with the corresponding Floquet multipliers by the relation 

where T = lix jio is the period of the multibreather. Due to the Hamiltonian character of the system there is a pair 
of exponents identically equal to zero. The non-zero characteristic exponents of the central oscillators correspond 
to the eigenvalues of the (2n x 2n) stability matrix E = VlD^H^^ , where VI is the matrix of the symplectic form 



-I 
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and I is the nx n identity matrix. The effective Hamiltonian H , as it has already been mentioned. 



in first order of approximation is given by — Hg + e{Hi). So, the stability matrix E, to leading order of 

approximation and by taking into consideration the form of H'^^ , becomes 



E 



A 


B 


C 


D 



eAi 




Co+£Ci 





— £ 



dcpidlj 



' dhdL 



—e- 



\ 



(5) 
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Since the only possible solutions are the ones with 0^ = 0, tt and we consider central sites oscillating with the same 
frequency uj, we get that Ai = Di ~ and so, the nonzero characteristic exponents are given to leading order of 
approximation by 

a±,; -±V£^ + 0(e3/2) (6) 

where xu *re the eigenvalues of the matrix Bi • Cq. Due to the form of the J I transformation the Cq matrix 
becomes (see [22| appendix B) 



Co 



dJ 



doj 
dJ 



( 2 

-1 





-1 



2 

-1 



2 / 
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So, (H]) becomes, up to leading order terms. 



± 



where Xzi are the eigenvalues of Z = Bi 
For systems of the form ^ we get 



Bi 



fi for i = j 
for i ^ j 



(7) 



with 



So, Z can be written as 



/ /i 
/2 



Z = Bi L = 



V 



/„_i 
/„ / 



m—1 



/ 2 -1 \ / 2/i 

-12-10 -/2 

0-12-1 

V -1 2 y V 



(8) 



-fi 

2/2 



— 1 2 — I — I 

-/„ 2/„ ; 



(9) 



Note that, for linear stability we require all the Floquet multipliers to lie on the unit circle, which is tantamount 
to all the characteristic exponents being purely imaginary. This depends on the sign of P = and the sign of Xz 
as it can be seen from ([T]). Finally, by using some counting theorems [l^l for ([SJ, we obtain: 

Theorem 1. f2^] In systems of the form {ip, if P = < the only configuration which leads to linearly stable 
multibreathers, for \e\ small enough, is the one with — tt Vz = 1 . . .n (anti-phase multibreather), while if P > the 
only linearly stable configuration, for \e\ small enough, is the one with (f>i — yi — 1 . . .n (in-phase multibreather). 
Moreover, for P < ( respectively, P > 0), for unstable configurations, their number of unstable eigenvalues will be 
precisely equal to the number of nearest neighbors which are in- ( respectively, in anti- ) phase between them. 

Remark 1: Note that, the form of the matrix Co is such due to the form of the J i— >■ / transformation and the 
fact that in the anti-continuum limit — lu and iOi = uj for i ^ 1 . . .n. So, it is independent of the range of the 
interaction between the oscillators of the chain and it will remain the same in what follows. On the other hand, the 
diagonal form of Bi will change if longer range interactions are added to the system. So, the theorem will no longer 
hold but the general methodology will still apply and the characteristic exponents of the multibreather will be given 
by ([7]). We will consider this case in what follows. 

Remark 2: In our previous works we used the term "out-of-phase" ioi (j> — tt configurations. This was because the 
only out-of-phase configuration was the </> = tt one. In the present work, since, as we will see in the next section, there 
are out-of phase configurations with (j) ^ n, we use the term "anti-phase" for the (p = n configuration. 



III. KLEIN-GORDON CHAIN WITH LONG RANGE INTERACTIONS 



The picture radically changes when the chain involves interactions with range longer than mere nearest neighbors. 
The range parameter r will be used to indicate the interaction length between the oscillators of the chain. So, for the 
classic nearest neighbor chain the range is r = 1 as shown in fig. [1] while for the next nearest neighbor (NNN) chain 
the range is r = 2 as illustrated in fig. [5] etc. The coupling force between the oscillators of the chain is linear and the 
coupling constants ei,i — 1 . . . r are not, in general, equal. 

The Hamiltonian of a ID KG chain with long range interactions is: 

oo 2 oo r 

E [y + ^(^«)] + 2 E E (^^ - ^■'+^-)' (10) 

2— — CXj i— — QO j — 1 

which leads to the equations of motion 

r 

Xi = -V'{xi) -f ^ - 2xi -I- Xi+j) 
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A. Persistence of multibreathers 



Let Ej = kjS, with ki = 1, then the Hamiltonian (ITOl) becomes 



oo 2 CO r 

Z — — oo 2— — oo J — 1 

Now, since the Hainihonian is written in the form Hq + e{Hi) the persistence conditions ([2]) can be used. If we 
consider again n + 1 "central" oscillators and Xi = X]m=o ^™ cos(rnz«i), we get for this case 



oo r n—j + 1 j — 1 

(^i)=-2EE E ^?„fc,cos(m^0,+,)- (12) 

m=lj"=l s=l i=0 

Note that, in the above we considered r ^ n since any interaction of oscillators with r > n does not affect the 
calculations, which are performed in the anti-continuum limit. So, if one considers r > n then for the calculations in 
this section it would be equivalent to the choice oi r ~ n. By differentiating Eq. [12] with respect to we get 

= ^ E E E E = 0' (13) 

or, by taking into consideration the definition of (|4]), 

Y,J^KM{J2<t>s+i) = 0, (14) 



where zi = max(l, i — p + 1) and Z2 



j for i + p — 1 ^ n 
n — p + 1 for i + p — 1 > n 

Eqs. [T2]and[T3l for[T4|) may be seem cumbersome to handle, they are much easier to use in the concrete examples 
that will follow in the next sections. 



B. Stability of multibreathers 

As we have already mentioned the characteristic exponents of the multibreather provided by the persistence condi- 
tions (fT3|l are given, to leading order of approximation, by ([7]), i.e. 



cr±t = ±\l -e^Xzi, i = l...n, 



where Xzi are the eigenvalues of Z with 



Z = Bi L 



d(f)id(j)j 



/ 2 -1 

-12-10 



v 



0-12-1 
0-12/ 



, i,j^l...n. 



(15) 



For linear stability we need all the characteristic exponents to be purely imaginary. So, if P = e^j < we need 



all the eigenvalues of Z to be negative, while if P = > we need all the eigenvalues of Z to be positive 
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FIG. 2: The r = 2 Klein- Gordon chain with next nearest neighbor interactions 

Without loss of generality we can consider i ^ j, since the Bi matrix is symmetric. Let d = j — i + 1, then the 



general form of 



IS 







if d> r 



d(j)id<f)j 



OG r Z2 



p-1 



m=lp=ds=zi 1=0 

or, by taking under consideration the definition of 



(16) 







p=ds=zi 1=0 



E E 

Now zi is given by zi — max(l, i — p + d), while Z2 is still given by Z2 



if d> r 
if d ^ r 



(17) 



i for i + p — 1 ^ n 
n — p + 1 for i + p — 1 > n 

In order to demonstrate the use of the results of this section, in what follows, we will examine some particular cases. 



IV. 3-SITE BREATHERS WITH r = 2 

A. The El — 62 = £ case 

1. Persistence of multibreathers 

The simplest case to check the effect of long range interactions is the one of 3 central oscilators (i.e. n = 2) and 
range r = 2. As we have already mentioned, any range r > 2 would not affect our calculations. First we will check 
the case fci = ^2 = 1 ^ £1 = £2 = £• 

In this case the Hamiltonian (fTTj) reads 

00 _^ 00 
H = Ho+eHi= + V{x^) + ^ ^ [{xi - x^+lf + {x^ - x^+2f] . 

2— — 00 Z— — 00 

Since we consider a 3-site breather. Hi becomes at the anti-continuum limit 

Hi = xf + xl + xl + {xi — X2)'^ + {xi — XzY + {x2 — XzY 

and 

00 

{Hi) = - 2 E {cos(m(/)i) + cos(to02) + cos[m{(j)i + 02)]} , 

771 — 1 



according also to (HU, for n = 2 and r = 2. The persistence conditions (|13l) become 
d{Hi 

- = U ^ 

in—l 



— = ^ mAl^ {sin(m(/)j) + sin[m(0i -f 02)1} = 0, for i = 1, 2 

"1 ■ f ^ 



(18) 



FIG. 3: Two snapshots of a 3-site (n = 2), anti-phase (01 = (^2 = t) muhibreather in a range r = 2 Klein-Gordon chain with 
e\ — ei — 0.02 and frequency a; = 2-k/7. See also videol in [ij. 

or, by taking into consideration the definition of M{(f) in (U), 

M{4>i) -f Af (01 + (1)2) = 0, for 1 = 1,2. (19) 
This equation, in addition to the standard solutions 

0j = 0,7r, 

provides also the solutions 

01 = 02 - 2^3, 01 = 02 = 47r/3. 

The multibreather solutions with 0^ 7^ are called phase-shift multibreathers or phase-shift breathers. The anti-phase 
and phase-shift configurations are depicted in figs. |3]and|4l For a better visualization one can also refer to videos 1 
and 3 in [1] (video2 shows an in-phase configuration). 

In order to produce these figures (and videos) we used the on-site potential V(x) = — 0.15 0.05— and 

initial conditions which correspond to motion with period T = 7 and frequency uj = 2'k/T — 27r/7 ~ 0.8976. The 
same potential is used for every numerical calculation throughout this work, although it is straightforward to apply 
the relevant notions to arbitrary potentials of the Klein-Gordon type. 
Remark 1: The persistence conditions ([19]) provide 2 equations. 



By substraction of equations ([20]) we get 



Af(0i)-t-M(0i+02)=O 
Af(02)+M(0i+02)=O. 



Af(0i) = M(02) (21) 

which has, besides the trivial solutions 0i = 0, tt, two other obvious solutions: 0i — 02 and 0i + 02 — tt for 
^ 0,; < 2tt. The last solution does not provide any new information because by substituting this into equations ([20)) 
we get Af(0i) = M{4>2) = 0, which , as it is shown in [20l|, only possesses the 0^ = 0, tt solutions. But the 0i = 02 = 
solution can reduce the two equations (|20p into equation 



A//(0) -f Af(20) = 0. (22) 



Remark 2: Our numerical computations strongly suggest that for all the phase-shift solutions it is 0i =02, yet a 
rigorous proof of this fact is still an open problem. So, equation <\22\ can be used in order to calculate all the solu- 
tions of the persistence conditions (|19p. except for the mixed one {0i = 0, 02 = tt} (or equivalently {0i = tt, 02 = 0}). 



Remark 3: In the case under consideration (n = 2, r = 2, fc^ = 1), all the available solutions correspond to 0i's 
which make each of the terms of the sum vanish in (jl8|) which obviously provides a zero total. 

Remark 4: The case under consideration is equivalent to the 3-site breathers on a hexagonal lattice which has 
already been studied in [l^Il^]. It can be effectively considered as a one-dimensional realization of such a lattice. In 
that context, the phase-shift multibreathers can be alternatively thought as "discrete vortices" , as they are solutions 
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FIG. 4: Three snapshots of a 3-site (n = 2), phase-shift [cj>\ = 4>2 ^ 0, vr) muhibreather in a range r = 2 Klein-Gordon chain 
with £1 = £2 = 0.02 and frequency u) = 27r/7. See also videoS in [![. 



which complete a phase rotation by 27r, as one traverses a discrete contour (which consists of the relevant triangle of 
sites). 

Remark 5: As an aside, it should be mentioned that an additional motivation for the consideration of such next- 
nearest neighbor interactions stems from the consideration of zigzag arrays, similar to the waveguide arrays proposed 
theoretically in the context of nonlinear optics (and hence in the realm of the DNLS equation) in [l^ • 
Remark 6: The stability of the above mentioned breathers will be discussed at the end of the next section as a 
special case of the more general unequal coupling one. 



B. The El 7^ £2 case 



1. Perststence of multibreathers 

Although the £i = £2 case is the easiest and allows us to perform some analytic calculations as well, the natural 
consideration for the case of next-nearest neighbors is the one with ei ^ £2- Intuitive physical considerations suggest 
to enforce £1 ^ £2 (considering coupling force decreasing with the distance between the oscillators) but there are 
configurations (like the zigzag one) which may also justify settings with £1 < £2 [l3l |. Let ki ~ 1 and k2 = k or, £1 = £ 
and £2 — ke. In this case, the Hamiltonian reads 

00 _^ 00 
H = Ho+eHi= -p^^+v{x^) + ^ ^ [{xi ~- x^-i)^ + k{x^ - . 

z— — 00 z— — 00 

Since we consider a 3-site breather (n = 2) we have only two independent (/)j's in the anti-continuum limit and by 
dm) we get 



^ 00 

{Hi) = - 2 X! {cos(m(/)i) -I- cos(m02) + A:cos[m(^i -I- (/12)]} • 

m— 1 

This leads to the persistence conditions: 



d{H_ 



— = 0^y^ mAl^{sm{m(f)^) + ksm[m{(j)i+ (t>2)]} = = M{(t)i) + kM {(1)1 + (j)2)^0 for i ^ 1,2. (23) 



m— 1 

Remark: By using the same arguments as in the previous section, if we consider (pi = (j)2, we get from 



^ mAl^ [am{m(f)) + k sin(2m(/))] = = A/ (0) + fcM(2(/)) = 0. (24) 



So, one could use p4|) instead of (|23p as the relevant persistence condition in order to calculate all the solutions of 

except of the mixed one {c/fi — 0, 4>2 = (or equivalently {0i = tt, 02 = 0}). 
In the fc = 1 (£1 = £2) case, one could make a choice of 4>i ^ 4>2 ^ 27r/3 or 47r/3 in order to have 
{sin(m(/>i) -|- sin[m(0i +(^2)]} = Vrn, so that the total sum in (jlSp would vanish also. This is not possible in 
the fc 7^ 1 case. So, one may be led to believe that this is an isolated solution and that possibly there are no other 
solutions than = and = tt in this case. However, it instead turns out that there can be other solutions also 
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FIG. 5: The bifurcation diagram for a 3-site (n — 2) breather in a r = 2 Klein-Gordon chain with k = e^jsi- A pitchfork 
bifurcation occurs for k — kcr ~ 0.48286. The curve is calculated using two methods. The first method is to use the full 
model, calculate the multibreather solutions with e = 0.001 and depict them as well as their stability. The second method is to 
solve numerically (|23[) and check when solutions with (pi ^ 0,tt appear. The curves produced with the two methods practically 
coincide (i.e., no difference is discernible at the scale of the plot). The various families that appear here are depicted in more 
detail in Fig. [G] 




(a) (b) (c) 

FIG. 6: The various families that constitute Fig. [5] are depicted. Portraits of the configuration of the central oscillators for 
the different multibreather families are shown as insets in the corresponding diagrams. The line type in the figures depend 



in the number of positive ( equivalent ly, the number of real eigenvalue pairs for P = ^-qJ < 0) the corresponding family 

possesses: no positive Xz corresponds to solid line, 1 positive corresponds to dashed line and 2 to dashed-dotted line. In (a) 
the in-phase configuration is depicted {(j>i = cf>2 = (or 27r)} which possesses 2 positive Xz. In (b) the mixed configuration is 
shown {(pi — 0, (j>2 = tt} which possesses 1 positive and 1 negative Xz. In (c) two families are shown. The first is the anti-phase 
one {(j)i = (j>2 = tt}. It has 2 negative Xz until k < kcr = 0.48286 while it has 1 positive and 1 negative Xz for k > kcr- At this 
point the anti-phase family bifurcates to provide the phase-shift configuration {(pi = 02 7^ 0,7r}. This family is represented by 
(pi ~ (j>2 ~ (2) or (pi — (p2 ~ (3) and has no positive (2 negative) Xz. AH of these together are depicted in Fig. [S] When two 
line segments coincide the more dense is shown. 

which can be calculated numerically for fc ^ 1. In fact, there is a critical value kcr — 0.48286 of k where a pitchfork 
bifurcation occurs (fig. [5]). For values k < kcr the only solutions Eq. ([23| [or ([M)) ] has are the trivial ones 0i = 0, tt. 
For k > kcr, i-G-, past the supercritical pitchfork bifurcation point, other solutions appear with (pi 7^ 0,7r (phase-shift 
breathers) as is shown in fig. [SJ 

The bifurcation curve has been calculated in two ways. Firstly by numerically modeling the full system and secondly 
by numerically solving the transcendental existence conditions using a small value of e = 0.001. The two curves 
practically coincide, which illustrates the remarkable accuracy of the theory in the vicinity of the anti-continuum 
limit. 

A phase-shift breather with k = 0.54 is depicted in fig. [T) For a better visualization of this breather one can also 
see video4 in [l|. 
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FIG. 7: Snapshots of a phase-shift 3-site breather for ei — 0.02 and £2 ~ feei = 0.56ei and frequency uj — 2tv/7. See also 
video4 in 

2. Stability of multibreathers 

By using the previously developed theory, we can calculate the characteristic exponents of the various configurations 
of 3-breathers in this lattice setting. The characteristic exponents of the specific solutions are given to first order of 
approximation by ^ as 

I K, 

where Xzi are the eigenvalues of the matrix Z defined in (jl5p . 

In the case under consideration of 3-site (n = 2) breathers with range r = 2 we have, also from (|17p . 

d'{Hi) ^ ( /(0i) + fc/(0i + (/.2) fc/(</.i + </.2) \ and L=f ^ 

^ kficj^i -f ^2) fih) + kf i^i +h) J [-1 2 J- 

So, we get from (fTS]) . 

2 ^ d^{H,) ^ ^ / 2/1 + fc/1+2 - /i \ 

d<l>id<l)j y kfi+2 - /2 2/2 + A:/i+2 y ' 

where the function /((/>) is defined as in and /1+2 = /(^ii -I- (/'2), while fi — f{4>i) for i — 1,2 . 

For linear stability it is required that all of the characteristic exponents be purely imaginary. So, the stability is 
determined by the sign of Xzi- 

In particular, we check the configurations that can appear in this case. 

• (pi =02- This is the general case and includes the in-phase {(pi — <p2 =0}, the out-of-phase {<pi — (f>2 — t^} 
and phase-shift configurations {(pi — (p2 ^ 0, tt}. The corresponding eigenvalues Xz are Xzi — 3/^ and 

Xz2 = f<t> + 2fc/20. 

• (pi =0,02 = TT. This is the only solution with (pi 7^ (p2. For this case it is Xzi2 — /o + (1 + fc)/7r i 
V/^-(l + A:)/o/. + /2(l-fc + fc2). 

Remark: We have that /o > as a direct consequence of the definition ([8]) of /(</>). On the other hand it is /^r < 0. 
This can be rigorously proven ([l^] Lemma 3) but it can also be intuitively understood by the definition (j8]) of j{(p') 
and the fact that the first term of the Fourier expansion of x{w) ([3]) is the dominant one. Using the same arguments 
we can conclude that |/o| > I/ttI- So, we can immediately conclude that in the in-phase {(pi — (p2 — ^{ configuration 
it is Xzi > 0, while in the mixed {<pi =0, (p2 = tt} configuration it is Xzi > 0, Xz2 < 0- 

On the other hand for the anti-phase {(pi — (p2 — t^} configuration, the formulas for the Xzi read 

Xzi = 3/^ and Xz2 = .U + 2k Jq. (25) 

The Xzi eigenvalue is always negative while the sign of the Xz2 depends on the value of k. This can provide us with 
a criterion about the value of k^r where the bifurcation occurs, since, at this point Xz2 changes sign. So, by (|25p we 

get Xz2 kcr = -TTT- 
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The values of /o and /tt depend on the particular on-site potential as well as on the frequency we examine, so, the 
value kcr is not fixed. But, if we consider breathers with relatively low amplitude, which amounts to the breather 
frequency uj being close to the phonon frequency ujp, a rough estimation of kcr can be made. In such a case, the 
nonlinear character of the system is not fully revealed yet which means that the Ai term in the development (jS]) is 
by far the most dominant one. This results to |/o| — I/ttI and consequently kcr — 0.5. 

In order to check our estimation, we perform some numerical calculations for the lattice with potential V{x) = 

— — '-'■^^"3 ^'^^^ which we use throughout this work, considering a motion with u = 2tt/7 ~ 0.8976. For this 

frequency, it is I/o/ZttI — 1, as can be seen in Fig. |8l so our estimation holds. In particular, it is /tt — —1.48658, 

/o ~ 1.53934 and kcr — ~777~ = 0.48286, which is precisely the value where the bifurcation occurs, while being very 

close also to the rough estimation (of 0.5) above. Note that, as it can be seen in Fig. [51 if we had chosen a smaller 
breather frequency uj, our estimation would be completely mistaken, since for small values of uj it is I/o/ZttI > 1- 




0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 8: Dependence of the \fo/ fw\ ratio with respect to the frequency cj of the breather. 



Remarks about figs. [5] and [6} In figs. [5] and [6] all the multibreather families that exist in the present configuration 
(n = 2, r = 2) are shown. The multibreather families correspond to solution families of Eqs. 1231 These families are 
categorized by the phase differences 4'i between the successive oscillators in the anticontinuous limit. The values of 

in the usual families = 0,7r) are constant with increasing fc, while in the phase-shift ((pi 7^ 0,7r) families their 
values change with respect to k. 

The various solution families are represented by various line (or curve) segments in the figures. The kind of the 

line depends on the number of positive Xz (i-e., of real eigenvalue pairs for P = ^'qJ < 0) that the corresponding 

solution has. So, for no Xz > we use a solid line, for one Xz > we use a dashed line while for two Xz > we use 
a dashed-dotted line. Since in Fig. [5] some of the families coincide, we separated the information in this figure into 
3 panels in Fig. [51 These 3 panels together compose Fig. [SI If the segments which represent two or more distinct 
families of solutions coincide, the more dense is shown in the figure. In order to facilitate the visualization of the 
various families, we added insets in Figs. [51 demonstrating the profiles of (and hence illustrating the phase difference 
between) the central oscillators in the anticontinuum limit. This has as a result only solid and dashed segments to 
appear in Fig. [5l The families that are depicted in the figure are: 

• {4>i — 4>2 — (or 27r)} (in-phase). This family is shown in Fig. [HJa) and possesses 2 positive Xz- 

• {4>i = or TT, (1)2 ^ tt} (mixed). This family is depicted in Fig.lHKb) and possesses 1 positive and 1 negative Xz- 

• {(1)1 = (/'2 = tt} (anti-phase). It is represented in Fig.[5Jc) by (j)i — (j)2 — (!)• It has no positive Xz until k < kcr 
while it has 1 positive and 1 negative Xz for k > kcr- At this point the 0i = 02 = tt family becomes subject to 
the bifurcation that gives rise to phase-shift multibreathers. 

• = 02 7^ 0, TT (phase-shift). This family is represented in Fig. M,c) by 0i = 02 = (2) or (f>i — (p2 ~ (3) and has 
no positive Xz- 

Since the stability of the multibreathers is also determined by the sign of P = £-^, the above are summarized, in 
terms of stability of the solutions, in Table [H 
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P 


k 


In-phase 


Out-of-phase 


Phase-shift 






<t>i = 4>2 = Q 


01 = 02 = TT 


01 = 02 7^ 0, TT 




P <0 




unstable 


stable 




Linear 


P < 


J\i J\i(2,T 


unstable 


unstable 


stable 


Stability 


P > 




stable 


unstable 






P > 




stable 


unstable 


unstable 



TABLE I: Stability of the various n = 2, r = 2, breather configurations depending on the values of P = ^ and k. With the 
dash we denote that this particular family does not exist for this range of values of k. 

Stability of the various 3-site breather configurations in the ei = 62 case: Using the above derived results 
we can conclude what it is already known from [2^^4|. i.e. for P < 0, as long as fc < kcr the only stable configuration 
is the anti-phase one, while for k > kcr the stable configuration is the phase-shift one (which corresponds in this case 
to the "vortex" configuration of [U the other hand for P > the only stable configuration is the in-phase 

one. 



V. 4-SITE BREATHERS WITH r = 2 

In the next configuration we will consider four central oscillators, in order to study larger configurations, but we 
will keep the range to r = 2 as a first step. 

A. Persistence of mult ibreat hers 

We will treat the two cases ei = £2 and ei 7^ £2 together, since the latter is a special case of the former with 
fci = ^2 — 1- Since we consider 4 (n = 3) central oscillators and range r = 2, (|12l) gives for fci = 1 and k2 = k, 

1 °° 

(Hi) = - 2 X/ {cos{m(f>i) + cos(rn02) + cos{m(/)3)+ 

m—l 

-ffccos[m((/)i -I- 02)] + fccos[TO(02 + 03)]} 
while, the corresponding persistence conditions ()13|) become 

M(0i) + fcAf(0i+02) = 

M(02)+fc[M(01+02)+M(02 + 03)] - 
M(03) + fcM(02 +03) = 

which have the trivial solutions (pi = 0, tt, as well as non trivial ones, as can be seen in fig. |9l 

By using the same arguments as in the previous section, which are also verified by our numerical investigation we 
have that for all the phase-shift breathers it is 0i = 03. At fc = fccr'' — 0.3219 the anti-phase {0i = 02 = 03 = tt} 
family becomes subject to a bifurcation that generates the phase-shift 4-site breathers. We should also note in passing 
(see details below) that, in addition to this supercritical pitchfork, the figure reveals also a sub-critical pitchfork 
bifurcation that terminates the two asymmetric branches upon their collision with the branch with the mixed family 
{01 = 03 = TT, 02 = 0} at ki? = 1.0736. 

B. Stability of multibreathers 

The stability of the existing multibreather solutions can be calculated by using the previously developed theory. 
Their corresponding characteristic exponents are given to first order of approximation by (I?]). In the case under 
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k 

FIG. 9: In this diagram all the existing families of 4-site (n = 3) breathers in a KG chain with range r = 2 are depicted. For a 
more detailed view of the particular families appearing in this diagram, refer also to fig. 1101 



consideration of 4-site (n = 3) breathers with range r — 2 we have from (|16p 



So, (US]) gives 



/l + fc/l+2 fc/l+2 

k fl+2 f2 + k fi+2 + k /2+3 k /2+3 
k /2+3 fs + k /2+3 




2/l + kfi+2 fc/i+2 - fi -kfi+2 

k{fl+2 - /2+3) - /2 2/2 + /s(/l+2 + /2+3) fc(/2+3 - /1+2) - /2 

-fc/2+3 fc/2+3 - /s 2/3 + fc/2+3 

For the general case 0i — 4>3, the eigenvalues of Z are 

X.i = 2(/i + fc/i+2) and X.2,3 = /i + /2 + fc/1+2 ± y//? + /I - 2A:/i/i+2 + ^2/2^3. 

The only configurations that are not included in the case above are the mixed ones {0i — (j)2 — 0, 4>z — tt} and 
{(j)i = 02 = TT, 03 = 0}, which both have 2 positive x^, independently of the value of k. 

Remark: The Xz2 eigenvalues of the anti-phase and mixed configurations can be used in order to calculate the 
values of fccr'' and fcir^ The Xz2 fo^' the anti-phase configuration is 



X.2 - 2/, + fc/o + ^2/2-2fc/,/o + fc2/2. 



1 y 1 

Since for fc = fccr'' it is X22 = we get k]^ — ^ ~ — . The last rough estimation can be performed only when 

3/0 3 

we consider breathers with frequency w close to the phonon frequency ujp (see also the discussion in the previous 
section), where |/o| ~ |/^|. In order to be more precise, for the potential and frequency used in the present work, we 
get k^r = 0.3219. 

On the other hand, the Xz2 eigenvalue for the mixed {<j)i = ^3 = tt, (j)2 ~ 0} configuration is 



(2) . . (2) 

Since for k = kcr it is Xz2 = we get kcr = — 



Xz2 = /o + (1 + k)U + v/(l - + /o 

/o 



work k, 



(2) 



2/7r + fo 



1 or, for the specific potential and frequency used in this 



1.0736. 
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(b) 



/\ 
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(c) 



0,,/7r 




0.0 0.2 0.4 



(d) 



FIG. 10: The various families that appear in Fig. are depicted. Portraits of the configuration of the central oscillators for 
the different multibreather families are shown as insets in the corresponding diagrams. The line (or curve) type in the figures 
depends on the number of positive Xz that the corresponding family possesses: no positive Xz corresponds to a solid line, 1 
positive Xz corresponds to a dashed line, 2 to a dash-dotted and 3 to a dotted line. In (a) two multibreather families are 
depicted. The first is the in-phase one with {0i = 02 = 03 = (or 2tt)} which possesses 3 positive Xz- The second one is the 
anti-phase one with {^i — 4>2 = 4>3 ~ It has no positive Xz until k < fccr^ = 0.3219 while it has 1 positive Xz for k > kil). 
At this point the anti-phase family gives rise, through a supercritical pitchfork, to the the phase-shift {(j>i = 02 7^ 0, vr} family. 
In (b) the mixed 1 configuration is shown {0i = 02 = or 27r, 03 — tv} or {0i = 02 = vr, 03 = or 27r}. These families possess 
2 positive and 1 negative Xz each. In (c) the mixed 2 family {0i = 03 = vr, 02 = or 2-k} is shown. It possesses 1 positive Xz 
for k < kiV = 1.0736 and no positive Xz for k > kir^ . At this point the mixed 2 family collides with the phase-shift one. In 
(d) the phase-shift {0i = 02 7^ 0, tt} family is shown and it is represented by 0i = 03 = (3) and 02 = (2) or 0i = 03 = (4) 
and 02 = (5). It exists for kcV < k < kiV ■ At fe = fe'r' it bifurcates from the anti-phase family, while at = fccr' it emerges 
with the mixed 2 family, and possesses no positive Xz- Fig. [9] contains all of the above families together, where, when two line 
segments coincide the more dense is shown. 



Remarks on the stability diagram of Figs. [9] and llOt In Figs. IH] and [TU] all the multibreather families they 
exist in the present configuration (n = 3, r = 2) are shown. The kind of the line (or curve) used for every segment 
depends on the number of positive Xz as follows: no positive Xz — > solid, 1 positive Xz — > dashed, 2 positive Xz — ^ 
dashed-dotted, 3 positive Xz — >■ dotted. All the families shown in FiglTU]are depicted together in FiglHl Since, if two 
segments coincide, the more dense is shown, in FigIS]we can see only solid and dashed segments. The families which 
are depicted in these two figures are the following: 



• {01 = 02 = '/'3 = (or 27r)} (in-phase). This family is shown in Fig llOf a) and has 3 positive Xz- 

• {01 — 4>2 ~ (j^s — ''^} (anti-phase). It is shown in Fig flUT a) and has no positive Xz for k < — 0.3219 while it 
possesses one positive Xz for k > kcl-^ . At this point the anti- phase family bifurcates to provide the phase-shift 



breather family. 



TT, 



0} (mixed 1). These families are depicted in Fig flUf b) and they 



• {01 = 02 = 0, 03 = tt} or {01 = (, 
both possess 2 positive Xz 

f 21 

L13 = n, 02 = (or 27r)} (mixed 2). It is shown in Fig iTUT c). For k < kcr it has 1 positive Xz, while for 

1.0736, this family collides with the phase-shift family. 



{01 = 03 

k > kcr it has no positive Xz- At k = kcr 
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• (pi — 03,02 7^ 0,7r (phase shift). This family is represented in FiglTUTd) by 0i = 03=(3) and 02 = (2) (or 
01 = 03=(4) and 02 = (5)) and has no positive Xz- It begins to exist at fc = kil) where it bifurcates from the 

(2) 

anti-phase family and cease to exist at k ^ kcr where it coUides with the mixed 1 family. 

Since the stability of the multibreathers is also determined by the sign of P = the above are summarized, in 
terms of the linear stability of the corresponding configurations, in Table HIl 





P 


k 


In-phase 


Out-of-phase 


Phase-shift 


Mixed 






01 = 02 = 03 = 


01 = 02 = 03 = 7! 


01 = 03, 02 7^ 0, TT 


01 = 03 = TT, 02 = 




P < 




unstable 


stable 




unstable 




P < 


1,(1) ^ h ^ 


unstable 


unstable 


stable 


unstable 


Linear 


P < 




unstable 


unstable 




stable 


Stability 


P > 


k kf^^rp 


stable 


unstable 




unstable 




P > 


f^CV — cv 


stable 


unstable 


unstable 


unstable 




P > 


k ^ kcr 


stable 


unstable 




unstable 



TABLE II: Stability of the various 4-site (n — 3) and range r = 2, breather configurations depending on the values of P = ^ 
and k. With the dash we denote that this particular family does not exist for this range of values of k. 



VI. 4-SITE BREATHERS WITH r = 3 

The natural way to extend our study in 4-site breathers is to consider range of interaction r — 3 (i.e., involving 
interactions with the 3 closest neighbors on each side of the chain) , in order for all the central oscillators to interact 
with each other. 



A. Persistence of multibreathers 

Bearing in mind that — eki and that ki = 1, (Hi) in this case becomes 

^ oo 

(Hi) = - 2 X/ {cos(m0i) -I- cos(m02) + cos(m03)-|- 

m— 1 

-Hfc2{cos[m(0i + 02)] + cos[r7i(02 + 03)]} + ks cos[m(0i + (j)2 + 03)]} , 
while the corresponding persistence conditions become 

M(0i) + fc2M(0i + 02) + A:3M(0i + 02 + 03) = 

M(02) + fc2 [M (01 + 02) + M(02 + 03)] + A:3M(0i 02 + 03) = 
M(03) + fc2M(02 -I- 03) + A:3Af(0i -f 02 + 03) = 0. 

For every ki, there exist the usual 0i = 0, tt solutions, as well as others as it can be seen in fig. 1111 By keeping k^ 
constant, we get various mono-parametric bifurcation diagrams with k2 as the parameter. Again, for all the phase 
shift configurations it is 0i =03. In fig. 1111 the bifurcation diagrams for two values of are depicted, fcs = 0.2 and 
ks = 0.4. We see that the value of fc2cJ.'' where the supercritical bifurcation occurs depends strongly on the value of 
/c3, while the value of ^2^^-' remains almost constant at k2^^^ — 1.075. The dependence of ^2^ with respect to k^ is 
shown in fig. [T^ 

Note that for k^ this case coincides with the r = 2 case (i.e., the latter is a special case example) and we 
retrieve the diagram of Fig. HI 
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B. Stability 

As it has already mentioned the stabihty of the multibreathers is determined by the sign of the eigenvalues Xzii of 
matrix Z By (HBl) we get 



d(l)id(j)j 
and since 



/ /l + + ^3/1+2+3 ^2/1+2 + ^3/1+2+3 ^3/1+2+3 

^2/1+2 + ^3/1+2+3 /2 + ^2(/l+2 + /2+3) + ^"3/1+2+3 ^2/2+3+^3/1+2+3 



V 



1+2+3 



^2/2+3 + ^'3/1+2+3 



2-10 
-1 2 -1 
0-12 



h + ^2/2+3 + fcs/i 



+2+3 



Z = 



we finally get 



( 2/1 + ^2/1+2 + fca/iH 



fc3/ 



1+2+3 



1+2 



+2+3 ^2/1+2 — /l 

^2/1+2 + ^3/1+2+3 - /2 - ^'2/2+3 2/2 + fc2(/l+2 + /2+3) ^2/2+3 + ^3/1+2+3 - /2 - ^2/1+2 
\ ^3/1+2+3 — ^2/2+3 ^2/2+3 ^ /3 2/3 + fc2/2+3 + ^3/1+2+3 



Its eigenvalues are, for 01 — ^3 and non-specific values of 0i, 

X.l = 2(/i + fc2/l+2) 



X^2,3 ^ h+ h + ^2/1+2 + kzh<t>i+'t>2 ± /f + /I - 2/C2/1/1+2 + fci/1+2 - 2^3/2/201+^2 + fc|/|0i+02- 

Although such analytical formulas exist and accurately predict the stability and bifurcations of the system, a clearer 
understanding emerges from the observation of the associated bifurcation diagrams (fig. 1111) . The diagrams present 
exactly the same solution families as in Fig. [HI but in this case the value of k2^^r ^ where the supercritical bifurcation 
occurs, is strongly affected by the value of fcs, while, the value of fc2cr\ where the subcritical bifurcation occur, s 
remains almost constant with k^'c) — 1-075 fFig. [T^. This indicates that, for the range of values of ^3 considered in 
this figure, the parametric interval of k2 (the strength of next-nearest neighbor interactions) over which phase-shift 
solutions exist narrows as k^ (the strength of interaction with the third-nearest-neighbors) is increased. 




0.0 0.2 0.4 6 0.8 1.0 1.2 1.4 2 04 06 8 1 1 2 14 

^2 ^2 



FIG. 11: In these panels the various 4-site (n — 3) multibreathers families in a chain with interaction range r = 3 are shown. 
In the diagram we have considered k2 variable, while fca = 0.2 in the left panel and fca — 0.4 in the right panel. The families 
depicted are qualitatively the same as the ones in the n = 3, r = 2 case. The only difference is that in the present case the 
value of A;2ir' where the supercritical bifurcation occurs depends strongly on fca while the value of ^2!?' where the subcritical 
bifurcation occurs remains almost fixed. 
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FIG. 12: The values of the two fc2cr, where the bifurcations occur, with respect to fca. Although fc2cr' depends strongly on fcs, 
A:2cr remains almost constant at k^cr — 1.075. 



VII. 4-SITE BREATHERS IN A 2D SQUARE LATTICE WITH r = 2 

We now turn our considerations to the case of a square lattice, as the one in Fig. [131 with nearest-neighbor 
interactions, not only with the horizontal and vertical neighbors, but with the diagonal as well. The latter interaction 
is assumed to have a strength £2 = kei (where ei = e will be taken to denote the coupling strength of adjacent nodes 
along the lattice axes). The Hamiltonian for this system is 



7 = J 
7 = 2 
7=1 
7 = 
7 = -l 
7 = -2 







-1 / = i=l 



FIG. 13: In the square lattice under consideration each oscillator is coupled with its neighbour not only in the horizontal and 
vertical directions but in the diagonal directions as well. So, every lattice site interacts with its 8 neighbours instead of the 4 of 
the classical KG square lattice configuration. In this setting, we consider 4-site breathers. Let the encircled oscillators in the 
figure be the central oscillators which are denoted as 0, 1, 2, 3 respectively 
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(26) 



t — — 00 J — — 00 



or 
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H = Ho+ sHi 



1 



E E 

■i— — oo j — — o 
oo oc 

f E E {( 



plj + V{xij) 



(27) 



1 — — 00 J — — QO 



We consider 4 "central" oscillators in the ant i- continuum limit and we denote them by 0, 1,2 ,3 as it is shown in fig |13l 
We have then 01 = wi — Wq, 4>2 — W2 — wi, (f>3 — — W2 and 4>i = wq — w^. But since, by construction, we have 
(^4 = 27r — — 02 — 03, we have finally only 3 independent 0i's. The (Hi) in this case is 

-j^ oo 

(Hi) = - 2 E "^'n {cos(m0i) + cos(m02) + cos(to03) + cos[m(0i + 02 + 03)] 

m— 1 

+fc{cos[m(0i + 02)] + cos[m(02 + 03)]}} 
and the corresponding persistence conditions are 

M(0i) + M(0i + 02 + 03) + kM (01 +02) =0 

M(02) + Af(0i +02 +03) +fc[M(0i +02) +M(02 +03)] = 
M(03) + A/(0i+ 02 + 03)]+/cAf(02 + 03) = 

This case coincides with the ID 4-sitc r ~ 3 chain case, with ~ ki — 1 and ^2 — k. Hence, the results (both 
for the persistence and the stability of the solutions) are a special case of the previous section. All the existing 
multibreather families of this configuration are depicted in fig. 1141 

In this diagram we can observe the appearance of {0i = 02 = 03 = ^/^l family, which is the vortex solution of the 
classical square Klein-Gordon lattice; see e.g., the relevant discussion in [ll|. In addition, there are several phase-shift 
families, the stability of which will be analyzed below. Interestingly, all these phase-shift breathers cease to exist at 
a critical value of fc = kcr — 1.03549 except of the vortex one. 
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FIG. 14: In this diagram, the various multibreather families of the square lattice configuration, are depicted. The value of k 
determines the strength of the coupling in the diagonal direction. For a more detailed description of the particular families one 
can also refer to Fig. 1151 

The mutibreather families that are supported by the present configuration are described in what follows. 
• {01 = 02 = 03 = 0} (in-phase) It is shown in FigfT^a). It possesses 3 positive Xz independently of the value of 

k2. 
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(c) (d) 

FIG. 15: The various multibreather families that consist Fig. [14] are depicted. Portraits of the configuration of the central 
oscillators for the different multibreather families are shown as insets in the corresponding diagrams. Note that, for better 
presentation, instead of showing the the 4 central oscillators in a square configuration, we show their ID equivalent. The 
line (or curve) type in the figures depends on the number of positive Xz the corresponding family possesses: no positive 
corresponds to solid line, 1 positive Xz corresponds to dashed line, 2 to dashed-dotted and 3 to dotted line. In (a) three families 
are shown. The first family is the in-phase one {4>\ = (j>2 = (f)^ = G ox 2-k}, which possesses 3 positive Xz. The second family 
is the vortex {i^i = (j!)2 = 03 = 7r/2 (or 37r/2)} family. It possesses no positive Xz. The third family is the anti-phase one 
{4>i = 02 = (/"s = tt}. This family possesses no positive Xz for k < kir^ = 0.96572, while for k > fcir' it acquires 2 positive 
Xz. At this point it bifurcates to provide the phase-shift 1 family. In this figure only the inset of the vortex family is present 
because all the others are the same as the ones presented for the equivalent ID configuration, in figs. 1101 In (b) the phase-shift 

1 {01 = 03 = TT, 02 = 0} family is shown, which exists for fcir' < k < fc'r' = 1.03549. It is depicted by 0i = 03 = (1) and 

02 = (2) or 01 = 03 = (1) and 02 — (3). It possesses no positive Xz for kil^ < k < fcir^ — 1.0344. At k — k^^ it collides with 
the phase-shift 2 family. For kcr < k < kcr the family possesses 1 positive Xz. At fc = fccr it collides with the mixed family. 
In (c) the mixed family {0i = 03 = tt, 02 = 0} is depicted. It has 2 positive Xz for k < kiV and 1 positive Xz for k2 > ki^K 
Finally in (d) two phase-shift families are depicted. The phase phase-shift 2a family, which is represented by 0i = 03 = (4) and 
02 =(8) which collides with the 0i = 03 = (5) and 02 =(6) for k2 = k^\ and the phase-shift 2b family which is represented by 
01 = 03 = (4) and 02 =(9) which collides with the 0i = 03 = (5) and 02 =(7) for = kiV ■ They both possess 1 positive Xz. 



• {4>i — 4>2 — <t>3 ~ ""/^l (vortex). It is shown in Fig ll5f a). It has no positive Xz independently of the value of 
and does not interact (i.e., collide) with any other family of solutions. 

• {01 = 02 = 03 = Ti"} (anti-phase). It is shown in Fig fTST a). It has no positive Xz for k < kil-^ = 0.96572, while 
for k > fccr'' it acquires 2 positive Xz- 

• {01 = 03 = TT, 02 = 0} (phase-shift 1). It is represented by in FigfTSTb) 0i = 03 = (1) and 02 = (2) or 
01 = 03 = (1) and 02 = (3). It exists for kil^ < k2 < kif' . It bifurcates from the anti-phase family at fc = fccr'' 
and possesses 3 negative Xz for kcl'' < k < k^ = 1.0344. At fc2 = fcir^ it collides with the phase-shift 2 family in 

(2) (3) f3) 

a bubcritical pitchfork bifurcation. So, for k^. < k < kcr it possesses 2 negative and a positive Xz- At k — kcr 
it collides with the mixed family. 
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• — 413 — TT, (j)2 — 0} (mixed). It is shown in Fig fTST c'). It has 2 positive Xz for k < kcr — 1.0355 and 1 
positive Xz for > kcr'. 

• (phase-shift 2a) It is represented in Fig fTST d) by (pi — 4>3 — (4) and 02 =(8) which coUides with the 0i = 03 = (5) 

(2) 

and 02 =(6) for k2 = kcr and have 1 positive Xz- 

• (phase-shift 2b) It is represented in Fig fTSl d) by 0i = 03 = (4) and 02 =(9) which coUides with the 0i = 03 = (5) 
and 02 =(7) for fc2 kcr and have 1 positive Xz- 

Notes on Figs. [14] and I15t 

1. At fc = kcr the 02 = (3) branch of the phase-shift 1 family coUides with the phase-shift 2a family, while the the 
02 = (2) branch of the phase-shift 1 family collides with the symmetric (and equivalent) of branches 02 — (7) 
and 02 = (9) of the phase-shift 2b family. 

2. At fc = the phase-shift 2 family approaches very much the vortex family, so it is plausible to expect that 
at that point they coincide. Yet, there is a small difference in the values of 0i of the two families due to the 
nonlinear character of the oscillators constituting the lattice, i.e., due to the existence of an infinity of terms in 
the development ([3]) . For a smaller frequency of oscillation, the nonlinear character of the oscillation becomes 
stronger, hence the terms after the first in ([3]) become larger and the two families are more clearly separated. 



VIII. DISCUSSION - COMPARISON WITH THE DNLS RESULTS 



It should be noted here that a number of results similar to the ones presented herein have been recently presented 
in the context of the DNLS equation e.g. in The setting of the DNLS essentially reflects a special case example 

of our Klein-Gordon calculation where instead of the existence and stability conditions reflecting a sum over all the 
harmonics due to the U(l) invariance of the underlying model, only the first harmonic is present. Nevertheless, the 
latter is sufficient to induce a number of the conclusions that we inferred herein. In particular, next-nearest neighbor 
interactions create phase-shift multibreathers (which were also parallelized to discrete vortex breathers in hexagonal 



lattices), as illustrated in [18[. As also shown in the same work, the long range interactions may drastically affect 



the stability properties of two-dimensional discrete vortices (in square lattices). On the other hand, the work of [3| 
provided a different analytical handle, via variational approximations, on the solutions that arise in settings with long 
range interactions. Furthermore, it was able to capture phenomena (both analytically and numerically) such as the 
supercritical or subcritical pitchfork bifurcations for such phase-shift multibreather solutions with NNN interactions. 
For instance, in the DNLS case the supecritical bifurcation leading to the emergence of such solutions would happen 
precisely at /c = 0.5 (due to the relevance of just the first harmonic) and not at fc = 0.48286 as obtained here in 
section IV B for the Klein-Gordon case. Nevertheless, the basic phenomenology remains intact. 



IX. CONCLUSIONS 



Classical Klein-Gordon chain with nearest neighbor interactions support multibreather solutions only with phase 
differences between successive oscillators of 0i = 0, tt. There, the stability scenaria are specific and well known. For a 
KG chain with P = < the anti-phase configuration is the only stable one, while for P = > the in-phase 
configuration is the only stable multibreather solution. 

On the other hand, in chains with long range interactions the picture is substantially different. First of all, in 
such chains, multibreathers with 0^ 7^ 0,7r (phase-shift multibreathers) can be supported in addition to the standard 
0i = 0, TT ones. The existence of phase-shift multibreathers as well as the specific 0i's of such profiles depend on the 
various coupling parameters Si within the chain. There are critical values of ki = Si/ei past which a bifurcation occurs 
(typically a supercritical pitchfork) and phase shift breathers begin to exist. Past this bifurcation point, the stability 
properties of the existing multibreathers are significantly modified, although this also depends on the particular (soft 
or hard) nature of the nonlinearity. As, however, additional parameters are tuned (e.g., higher ranges of neighbor 
interactions), it is also possible for such phase-shift solutions to terminate in subcritical pitchfork bifurcations. 

These results are not unique to the realm of one-dimensional lattices with higher range of interactions. They can 
also be developed for two-dimensional square lattices in which case they may lead to bifurcations or terminations of 
the families of discrete vortices which arise therein. Such vortices sustained by the two-dimensional analogs of the 
lattice can be of either a symmetric or asymmetric type. In particular, in the case considered herein, the presence of 
diagonal coupling within the square was critical in inducing the emergence of asymmetric such patterns. 
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This study opens a number of a directions for further investigation. Firstly, it would be very relevant to examine 
particular functional forms of the decay of the long range interactions (e.g., exponentially or polynomially decaying 
ones) to identify whether any systematic conclusions can be derived on the basis of such decay laws. Secondly, it would 
also be very interesting to examine the interplay of the geometry of higher dimensional lattices (and the interactions 
that they present) with the strength of the long range interactions that can be considered therein and to try to derive 
some general conclusions about the possible stable/unstable discrete soliton and discrete vortex solutions. Finally, an 
important and immediate direction that can be followed with the results of the paper in hand could be the effect of 
long-range interaction in phase-shift phonobreathers, whose stability for nearest-neighbor interaction was considered 
in [10| . A physical application of relevance and worthwhile of further investigation concerns the biological models for 
DNA or protein alpha- helices where dipole long-range interactions are relevant. 
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